Low self-afRne exponents of fracture surfaces of glass ceramics 
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The geometry of post mortem rough fracture surfaces of porous glass ceramics made of sintered 
glass beads is shown experimentally to be self-afBne with an exponent ( = 0.40 ± 0.04 remarkably 
lower than the 'universal' value = 0.8 frequently measured for many materials. This low value of 
( is similar to that found for sandstone samples of similar micro structure and is also practically 
independent on the porosity (j> in the range investigated (3 % < < 26 %) as well as on the bead 
diameter d and of the crack growth velocity. In contrast, the roughness amplitude normalized by d 
increases linearly with (j) while it is still independent, within experimental error, of d and of the crack 
propagation velocity. An interpretation of this variation is suggested in terms of a transition from 
transgranular to intergranular fracture propagation with no influence, however, on the exponent 

PACS numbers: 62.20.Mk, 46.50.+a, 68.35.Ct 



Since the pioneering work of Mandelbrot [l| , statistical 
mechanics has been used to describe fracture processes. 
A special attention has been brought to the statistical 
properties of the roughness of fracture surfaces which of- 
ten display a scaling behavior referred to as self- affinity. 
This means that they are statistically invariant under the 
transformation : r — > Ar and h — > X'^h; (r characterizes a 
location in the mean plane of the fracture, h is the sur- 
face elevation and C is the self-affine exponent Bi). Due 
to this invariance, the height-height correlation function 
Ah{Ar) =< {h{r + Ar) - h{r)f scales like : 



Ah 
Ah{d) 
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Ar is a vector of constant orientation and of modulus Ar; 
d will be taken equal to the characteristic length scale of 
the microstructure of the material (the mean grain size 
for instance). As a consequence, the amplitude of the 
roughness may be characterized by the value of Ah{d) 
(amplitudes corresponding to other Ar values can then 
be obtained from Eq. 

Experimental studies of various fractured 31? materials 
provided similar values of roughness exponents C, ~ 0.8 
(see Ref. Q for a review). These results supported the 
idea of a universal roughness exponent C independent on 
the material. However, sandstone fracture surfaces were 
found to be self-affine with a lower exponent close to 0.5 
0. Explaining these discrepancies is important to un- 
derstand the physical origin of their self-affine geometry. 
The objective of the present paper is to analyze these 
effects on an artificial material comparable to sandstone 
but with a microstructure that can be varied in a con- 
trolled manner. 

Glass ceramics made of sintered glass beads have been 



selected for that purpose because, like sandstone, they 
are made of cemented grains. However, in contrast with 
natural rocks, both the characteristic size and the co- 
hesion of the grains can be adjusted by modifying re- 
spectively the bead diameter and the sample porosity <f>. 
The exponent C for these materials will be shown to be 
much lower than the "universal value" 0.8; in addition, 
it depends surprisingly little on the microstructure of the 
material, although the amplitude of the roughness varies 
significantly with the porosity. The influence of velocity 
on roughness reported by some authors Qwill also be 
evaluated by comparing two types of mode I fracturing 
processes with different crack propagation velocities. 

The samples are prepared by heating a mold filled with 
a random packing of beads at 700° C during 20 to 200 
minutes. Depending on the duration of the annealing, the 
degree of sintering varies and porosities cf) ranging from 
3 % to 26 % are obtained. Either of two sets of beads 
with diameters in the range 104 — 128 /xm and 50 — 65 ^m 
are used. In the following, the characteristic microscopic 
length scale d introduced in Eq. is taken equal to 
the mean bead diameter. The open porosity is measured 
by saturating the sample with water and the profile of 
the total porosity along the sample is measured by 7- 
ray absorption. Porosity variations along the samples are 
of the order of 1 % so that (/) is considered as constant 
and equal to this mean value within ±1%. The process 
provides sintered glass cylinders with a 130 mm height 
and a 40 mm radius from which the samples used in the 
fracture tests are cut out. 

Two mode / fracture tests have been performed. Fast 
fracture propagations are obtained by means of modified 
Brazilian fracture tests using toroidal samples (inside and 
outside radii : 10 mm and 40 mm respectively; height par- 
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allel to the axis : 30 mm). A uniaxial radial compression 
is applied to the torus and two symmetrical dynamical 
cracks propagate from the central hole toward the outer 
zones where the compressive forces are applied. 

Quasistatic fracture propagations are realized on sam- 
ples with a modified Tapered Double Cantilever Beam 
(TDCB) geometry : Their width and length are respec- 
tively 20 mm (perpendicular to the crack propagation) 
and 60 mm (parallel to it). The fracture is initiated from 
a straight notch (thickness 1 mm) by applying on both 
sides a uniaxial tension with a constant opening rate : 
The tapered shape of these specimens allows us to obtain 
a stable mode I crack growth (see Ref. 0). The crack 
propagation velocity Vf, determined from the variations 
of the electrical resistance of a thin gold layer deposited 
on the side of the sample, remains close to 3 mm.s~^. 
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FIG. 1: Height profiles extracted from fracture surfaces of 
three sintered glasses with three different porosities but with 
beads with the same diameter 104 — 128 /im. 

Fig. ^displays fracture surface profiles of three samples 
of different porosities. The roughness amplitude increases 
by almost two orders of magnitude from ~ 1 /.tm for the 
densest {(j) — 3 %) to 100 ^m for the most porous 
sample {(f> — 26%) : Different profilometers are there- 
fore needed to scan all the samples. For surfaces with a 
very low spatial roughness (and therefore porosity) we 
use an interferometric optical profilometer ("'"'^Talysurf 
CCI 6000) with a vertical resolution better than 0.1 nm 
and a lateral resolution ~ 1 /xm. The local slope of the 
surface must always be small to apply this technique 
so that, for ^ > 7 %, a mechanical stylus profilometer 
(™Talysurf Intra) with a ~ 10 nm vertical and a ~ 2 /im 
lateral resolutions had to be used. For lower porosities, 
both profilometers are usable which allowed to check the 
consistency of the two measurements. For both profilome- 
ters the maps contain 1024 x 1024 points and the fields of 
view are respectively 3x3 and 6x6 mm. In the second 
technique, the stylus is always in contact with the surface 
and, for porosities higher than 18%, it often gets jammed 
into the deepest asperities of the surface. One uses then 



a point by point mechanical profilometer : A sensor tip is 
lowered down until it touches the surface in order to mea- 
sure its height; the tip is then raised by 200 ^m before 
getting moved laterally by 25 ^m to the next measure- 
ment point. The vertical and the lateral resolutions are 
respectively ~ 3 /.tm and ~ 10 /xm and the typical field of 
view is 6 X 6 mm. The consistency of the measurements 
was verified by comparing profiles provided by the two 
mechanical systems for (j> = 18 %. 

These measurements provide surface elevation profiles 
as shown in Fig. ^ In the following, we shall only ana- 
lyze profiles parallel to the z-axis {i.e. to the crack front) 
and located far enough from the initiation region so that 
the roughness properties are statistically stationary (this 
latter point is discussed in more detail below.) 
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FIG. 2: Log-log representation of Aft (averaged over profiles 
at different x values) as a function of Az for surface roughness 
profiles of several sintered glass samples. - Samples fractured 
using the TDCB test with porosities = 7 % (x) , 15 % (□) 
and 26 % (-I-) - Samples fractured using the Brazilian test 
with porosities ^ = 3 % (o) and 25 % (o) - Range of bead 
diameters used to realize samples : 104 — 128 fim - Straight 
lines : Linear fits with slope ^ (Tab.HJ. 



Quantitatively, these profiles will be characterized by 
their ID height-height correlation function, A/i(Az) =< 
{h{z + Az) — h{z))'^ : For surfaces with a self-affine 
roughness, the power law relationship should be ver- 
ified. This relation is tested in Fig. |21 which displays in 
a log-log scale the variations of Ah{Az) as a function of 
Az for several samples : These have different porosities 
in the range : 3 %<</>< 26 % and have been fractured 
both in a quasistatic mode and by the modified Brazil- 
ian test. The curve corresponds to an average of Ah{Az) 
over profiles corresponding to different distances x to the 
initiation and located in the region where the statistics 
of the roughness is stationary. 

Let us examine, for instance, the lower curve cor- 
responding to a 3 % porosity sample fractured using 
the modified Brazilian test procedure. The variation can 
clearly not be fitted by a single power law over the full 
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range of values investigated : Three domains of varia- 
tion are visible and correspond to exponents respectively 
equal to 1, 0.36 and (straight lines). The last value in- 
dicates that the surface appears as a plane at large scales 
and the first one corresponds to an Euclidian geometry 
of individual grains. The surface profile is therefore self- 
afHne (here with an exponent C = 0.36) only in the inter- 
mediate domain between two limiting length scales : In 
the log-log plot of Fig. 12 they correspond to the inter- 
sections between the straight lines fitted in the different 
domains. The lower boundary is of the order of the bead 
radius d/2 and the upper one will be referred to as Lc- 

This result is generalized by comparing the different 
curves in Fig. |21 corresponding to samples of different 
porosities and fractured both in the quasistatic and fast 
propagation modes. All the curves have the same global 
shape and, in log-log coordinates, the slopes are nearly 
the same in the intermediate domain : This shows that 
the roughness exponent is very similar in all cases while 
the vertical shift between the curves reflects different 
roughness amplitudes. 



similar limitation of the self-affine range of length scales 
by finite size effects has recently been reported 01 in ex- 
periments on mortar samples of different widths. 
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TABLE I: Physical and statistical characteristics of the sam- 
ples - (0) sample porosity; glass beads diameter range (B.D.I) 
104 — 128 /im, (B.D.2) 50 — 65 /im; crack propagation mode 
{Dyn.) dynamic, {Q.S.) quasi static; surface measurement 
technique (1) interferometric, (2) stylus profilometer (3) point 
by point - Statistical characteristic parameters of ID profiles 
normal to crack propagation (^ , ^ft) self-affine exponent val- 
ues obtained respectively from the variation of Ah with Az 
and from the Fourier power spectrum; {Ah{d)) roughness am- 
plitude in /xm; (Lc) upper boundary of self-affine domain in 
mm. 



The numerical values of the parameters characterizing 
all these curves are listed in Table for all samples in- 
vestigated in the present work. Table ^ confirms that ( 
has a very similar value C = 0.40 ± 0.04 for all samples 
independent of the bead size, of the porosity and of the 
crack propagation velocity: This common value is much 
lower that the value 0.8 reported for many materials and 
closer to the value 0.5 obtained for sandstone. This result 
is robust with respect to the method used to determine 
^ as shown by the comparison with the values (^ft in 
Table n obtained from the analysis of power spectra 0. 
The upper boundary Lc of the self-affine domain is of the 
order of 1 — 2 mm i.e. about a tenth of the sample width 
and seems to decrease for faster fracture propagations. A 
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FIG. 3: Variations of the roughness characteristic exponent 
and amplitude Ah{d) with the distance x to the initiation 
notch for the same sample (0 = 3%) as in FigOl 

In the above interpretations, it is assumed that, after 
a large enough propagation distance from the initiation, 
the statistical properties of the roughness become con- 
stant. In order to verify this assumption, the transient 
regime already reported by other authors [1,0] is studied 
in Fig. 13 In this figure, data points correspond to sin- 
gle profiles parallel to z at a given distance x from the 
initial notch. The roughness of each profiles is character- 
ized by its amplitude Ah{d, x) and by the exponent C,{x). 
For X < 1 mm, a transient regime in which Ah{d, x) de- 
creases with the distance x while Ci^) remains constant 
is indeed observed. At larger distances, ({x) and Ah{d, x) 
merely fluctuate around an average value. For the other 
samples investigated, the distance has the same order of 
magnitude. These results allow to restrict (as above) the 
statistical analysis to distances > 1 mm where the rough- 
ness is statistically stationary : This justifies our use in 
Fig. 121 of curves obtained from averages over several pro- 
files corresponding to different x values. 

The curves of Fig. |21 have confirmed the extreme ro- 
bustness of the exponent C which remains constant even 
in the transient regime. In contrast, the roughness am- 
plitude Ah{d) appears to be more sensitive, not only to 
X, but also to the characteristic parameters of the sample 
like its porosity as seen in Fig. ^ This variation is an- 
alyzed quantitatively in Fig. 0] in which the normalized 
variable Ah{d)/d is plotted as a function of the porosity 
(j>. For low porosities, Ah{d)/d decreases to zero while it 
is close to 0.5 for the most porous sample and a good 
collapse of all data on a linear variation is observed. This 
suggests that the velocity Vf of the crack propagation has 
a weak influence on the amplitude Ah{d) of the rough- 
ness and that Ah{d) is roughly proportional both to the 
bead size d and to 0. 

These results and, particularly, the increase of Ah{d) /d 



4 



2, 

< 0.2 
0.1 




10 20 30 

(%) 

FIG. 4: Variation of Ah{d)/d as a function of porosity <j). - 
Range of bead diameters used to realize samples -104— 128 /im 
(A) : Dynamic fracture (□) : Quasistatic fracture - 50 — 65 /im 
(0) dynamic fracture. The slope of the dashed line is 1.7. 



als using a 2D analysis based on the computation of a 
2D correlation function : Applying this latter technique 
to the present samples might help confirm the presence 
and amplitude of the anisotropy. These results may be 
compared to the theoretical predictions of brittle fracture 
models describing crack growth as the propagation of an 
elastic line in random heterogeneous media lisl , 
leading to C = 0.39 in the 2D case 13 El and the 3D 
cas e Il4l . Ductile fracture models such as those developed 
in |l7| and leading to higher roughness exponent values 
fail to reproduce the observations reported here because 
damage processes are not expected during the failure of 
glass ceramics. Further studies are currently in progress 
to confirm that these models are applicable to such ma- 
terials as glass ceramics or sandstone displaying brittle 
fracture and whether the higher value ^ — 0.8 measured 
for many others may be associated to ductile fracture. 



with (j) may be related to phenomena at the scale of a 
bead diameter. For high porosity samples, cracks propa- 
gate by breaking cemented necks binding two beads : The 
difference in height between neighboring beads is then of 
the order of their radius and Ah{d)/d ~ 0.5. For low 
porosity samples, neighboring beads are more strongly 
welded to each other and cracks propagate through the 
beads : The deflections of the surface are weaker com- 
pared to the bead radius and Ah{d)/d is lower. This is 
consistent with the increase of Ah{d)/d with </> in Fig. 
01 These differences in the propagation of the cracks are 
confirmed by scanning electron microscope (SEM) im- 
ages of fractured samples : These display a transition 
from transgranular to intergranular propagation as the 
porosity increases. It is remarkable that this transition 
has no influence on the characteristic exponent C- 

To conclude, the fracture surfaces of sintered glass sam- 
ples analyzed in the present work all display a self-affine 
domain of length scales characterized by an exponent 
C = 0.40 ± 0.04 : This value of C has been found to be 
independent of the porosity of the fracture propaga- 
tion velocity Vf and of the diameter d of the beads. C is 
significantly lower than for many materials {( ~ 0.8) but 
close to the value found for sandstone {( ~ 0.5) which has 
a similar structure. This nearly constant value of C. con- 
trasts with the linear increase with (f> of the normalized 
roughness amplitude Ah(d)/d (although it is indepen- 
dent both of Vf and d). This increase of the amplitude 
likely reflects a progressive transition from transgranular 
to intergranular fracture propagation also demonstrated 
by SEM observations : A surprising feature is the fact 
that this transition docs not influence Another impor- 
tant issue is the anisotropy of the surface roughness : Pre- 
liminary measurements of profiles parallel to the fracture 
propagation have given a value C. — 0.5 larger than that 
reported above for profiles parallel to the front. Similar 
differences have been reported 0, 0| on other materi- 
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